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The interaction of two spatially separated impurity atoms through phonon exchange in a Bose- 
Einstein condensate is studied within a Bogoliubov approach. The impurity atoms are held by deep 
and narrow trap potentials and experience level shifts which consist of a mean-field part and vacuum 
contributions from the Bogoliubov-phonons. In addition there is a conditional energy shift resulting 
from the exchange of phonons between the impurity atoms. 
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I. INTRODUCTION 



The ability to engineer the collisional interaction of 
ultra-cold individual atoms or ions as well as degenerate 
ensembles of atoms, such as Bosc-Einstein condensates 
(BECs) Q has dramatically improved in the last cou- 
ple of years by the development of quantum-optical tools 
such as single-atom micro-traps 0, La, , optical lattices 
0- 03- ■ atom-chips @ and others. Controlled collisional 
interactions of individual atoms are of fundamental in- 
terest but also have important potential applications in 
quantum information processing Recently the cou- 
pling of single-atom quantum dots to Bose-Einstein con- 
densates was studied in [ljj and the use of an impurity 
atom in a 1-dimensional optical lattice as an atom tran- 
sistor was proposed 0. We here study the mutual in- 
teraction between two separated, well localized impurity 
atoms through the exchange of Bogoliubov phonons in 
a BEC at zero temperature. When the impurity atoms 
undergo a state-dependent scattering with the conden- 
sate atoms, in addition to mean-field level shifts and lev- 
els shifts from the interaction with the vaccum fluctua- 
tions of the Bogoliubov phonons a conditional level shift 
emerges which results from phonon exchange between the 
impurities. This conditional shift is calculated and its 
dependence on trap geometry, impurity separation and 
the strength of the interactions within the condensate is 
studied. 



In section II we derive an effective coarse-grained in- 
teraction hamiltonian for the impurity atoms and relate 
the level shifts to correlation functions of quasi-particlc 
excitations. These will then be calculated within a Bo- 
goliubov approximation for a condensate in a box poten- 
tial in section III. It is shown that the coupling between 
the impurity atoms is strongest for a highly asymmet- 
ric geometry. For this reason we consider in section IV 
a quasi-one dimensional condensate. A simple analytic 
expression for the level shift is derived using a Thomas- 
Fermi approximation. 



II. EFFECTIVE INTERACTION OF IMPURITY 
ATOMS IN A BEC 

We here consider a Bose-Einstein condensate at T = 
with impurity atoms at fixed locations, which can be 
realized e.g. by tightly confining trap potentials as shown 
in Fig^ The traps are seperated such that any direct 
interaction of the atoms can be excluded. The atoms are 
assumed to have two relevant internal states |0) and |1) 
and shall undergo s-wave scattering interactions with the 
atoms of the BEC if they are in state 1 1) . If the traps are 
sufficiently deep, the atoms will stay in the corresponding 
ground state aV In this case the interaction hamiltonian 
of the condensate and the impurities has the form 
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dr |0o (r - ri)| ^ (r)tp(r) 
dr|0o(r-r 2 )| 2 ^ t (r)^(r; 



(1) 



where \a, (3) denotes the a-th internal state of the first 
and the /3-th internal state of the second impurity atom. 
The coupling to the condensate is described by the state 
dependent coupling constant n a with kq — and K\ = 
k. The condensate wave- function is denoted by tp. The 
ground state function of the impurities is given by 



0o (r) = 



: exp 



3,3 



2zl 



(2) 



with zo = yJh/msujQi ms being the mass of the impurity 
atoms and coq the frequency of the confining traps. 

In order to derive an effective Hamiltonian for the two 
impurity atoms it is convenient to first separate the in- 
teraction Q into a mean-field and a fluctuation part 



where 



|l)u(l 



^(^ + 11)22(11-^) 

a, (3 



(3) 



Q(t)= /drltoCr-rOWr,*)^,*), (4) 
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FIG. 1: Impurity atoms held by tight confining potentials in 
a Bose-Einstein condensate. When in internal state |1) the 
atoms undergo s-wave scattering interactions with the con- 
densate. 



and 



Bt(t) = d(t) - (C,(t) 



(5) 



The terms in the first line of eq. result in a mean- field 
level shift of the internal state 1 1) . They are of no interest 
in the present discussion and will be absorbed in the free 
Hamiltonian of the impurity atoms. 

We proceed by deriving an equation of motion for the 
statistical operator of the impurities interacting with the 
BEC. Within the usual Born approximation and as out- 
lined in Appendix A one finds 
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&0oi,oo - -T£a / dt'g 0l!00 (t')(B 2 (t)B 2 (t' 



dt 



' enM*){(Bi(t)BiW)) ( 8 ) 
Bi(t)B 2 (t')) + terms with 1 <-► 2 } . 



where the tilde denotes quantities in the interaction pic- 
ture and the matrix elements of the statistical operator 
arc denoted by Qap^s = (af3,t\ Q\j8,t). The correlations 

^BiBi'^ are calculated using the standard Bogoliubov 

approach, i.e. by setting 



(9) 



with tpQ being the solution of the Gross-Pitaevskii equa- 
tion and £ a small operator- valued correction and neglect- 
ing higher-order terms in £ (see Appendix B). Within 
the Bogoliubov approach we disregard terms of the order 

O (£ 4 ) in (b.bS and find 



(10) 



Bi(t)B v (t')\ = J2'^ iEi(t - t ' ) S j (l,l'). 



The Ej's are the Bogoliubov energies and 
S j (l,l') = Jdr |0o(r-r ; )| 2 Vo(r)K(r)~^(r)) 
x / dr' \M*' - r,)f Vo(r')K( r ') - • 



(11) 



The functions Uj and Vj are the solutions of the 
Bogoliubov-de Gennes equations (cf. Appendix B) and 
the prime at the sum indicates that the ground state is 
excluded. 

A calculation of the correlation functions shows that 
the often used Markov approximation cannot straight- 
forwardly applied to eqs.JUHJ). Instead we first use a 
Laplace transformation. Setting to — we find 



£ [Qa/3ns(t)] (p) 



0a/3,7<5(O) 



P + ■tfpMaprflp) 



(12) 



with 

r 1 v p+k E i 




P + 7L E i P-i E 3 J (13) 



P 



7i E 3 



p 



fi E 3 



K 2 - KpKS K 2 S - KpK S \ 

p+i E i p-k E i) 



In general, the Laplace transformation 11211 cannot be 
inverted analytically. However, if we are interested only 
in a coarse-grained time evolution, it is possible to ne- 
glect the p-dependence of M a p yl s, which amounts to 
M a f) n $(p) — ► M a p n g(0). In the coarse-grained picture 
the interaction of the impurity atoms with the conden- 
sate simply results into level shifts, i.e. 



The corresponding frequencies read 



(14) 



[Sj(l, 2) + Sj(2, 1)] (k 7 k 5 - n a Kp)} 



(15) 



This corresponds to an effective - coarse-grained - Hamil- 
tonian 



H cB = |10) (10| ftw 10 ,oo + |01) (01|nwoi,oo 
+ 111) (111 Woo. 



(16) 



The energy scheme of this Hamiltonian is shown in fig- 
ure |21 One recognizes from (|l(jf) for symmetric impurity 
locations a level shift 



^10,00 — ^01.00 



,.2 



Ah^ J 

j 



Sj{l,l) < 



(17) 
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of each impurity atom independent of the presence of 
the other. This level shift is due to the interaction 
with vacuum fluctuations of the Bogohubov quasi par- 
ticles (phonons). In addition there is a conditional level 
shift due to the exchange of Bogohubov quasi particles 
(phonons) between the two impurities: 



A 



^11,00 — w 10,00 

- 2 — < 1 
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FIG. 2: Energy scheme of the effective Hamiltonian for sym- 
metric arrangement of impurity atoms. Here a negative sign 
of A was assumed although positive values are possible. 



It should be noted that the coarse-graining approxi- 
mation is consitent with the collective level shift only if 



A <C - min' (Ej 



(19) 



where the prime indicates that the ground state is ex- 
cluded. In the following we will explicitly calculate the 
levels shifts for a homogeneous condensate, for an ideal 
condensate in a harmonic trap and a weakly interacting 
condensate in a trap in the Thomas- Fermi limit. 



III. HOMOGENEOUS CONDENSATE 

In this section we calculate the energy shifts 5 and 
A for the case of an interacting, homogeneous conden- 
sate with periodic boundary conditions of spatial peri- 
odicity L x , L y , L z , respectively. The solutions of the 
Bogoliubov-dc Gcnnes equations arc then given by plane 
waves w k = + / k )/2 and v k = (f£ - / k )/2 with 




ikr 



(20) 



The wave vectors k have to be chosen in such a way that 
they fulfill the periodic boundary conditions. Since the 



itk and i>k's have to be orthogonal to the ground state 
the case k = is excluded. The Bogohubov energies are 
given by 



>gn ) 



(21) 



with e k = h 2 k 2 /2mQ. By extending the integral in equa- 
tion 111|) over the whole M 3 , which is possible due to the 
effective cut-off provided by the impurity state wavefunc- 
tions <t>o, one can easily calculate the correlation func- 
tions: 

'Mt)B V {t')) = E'4 ex p {-{w-*) 

k 



x exp ik • (r, - v v ) - 



"0 K 



(22) 

Here No denotes the number of atoms in the condensate 
and V = L x L y L z . With this one finds 



K 2 Nn 



8 = T£ exp 



AhV 2 ' 

K 2 N - 

2hV 2 ■ 



E |^ cos ( k ■ Ar ) exp ( ~ 



(23) 



,(24) 



where Ar = V\ — V2- The sum over the Bogohubov quasi 
momenta converges due to the exponential term which 
effectively cuts off momenta with k 3> 1/zo- Obviously 
26 = A for ri = r2 and |A| < 2\S\. The conditional en- 
ergy shift A is shown in figurc|3J For very small distances 
of the impurities A is negativ and its absolute value ap- 
proaches its maximum, i.e. that of 25. For increasing 
distance the value of A increases monotonously and even- 
tually changes its sign. The monotonous increase would 
correspond to an attractive force between the impurity 
atoms if they could move freely. One recognizes, that 
for larger values of the dimcnsionlcss interaction param- 
eter K ~ g the energy shift decreases and the spatial 
dependence becomes less pronounced. This can be ex- 
plained by the increasing self-energy of the Bogohubov 
excitations. 

It is also very instructive to consider the dependence of 
the conditional level shift A on the condensate geometry, 
i.e. on the ratio £ ra d/^z, where X ra( j = L x = L y . This is 
illustrated in figure 0] One recognizes that the absolute 
value of the energy shift increases as the ratio L ra ^/L z 
decreases. Thus the energy shift is largest for a highly 
non-symmetric geometry of the BEC. The strongest ef- 
fect is thus to be expected in a quasi one-dimensional 
condensate. For this reason we will investigate in the fol- 
lowing section the energy shift in the case of a BEC in a 
harmonic trap only for a one-dimensional condensate. 



IV. 1-D CONDENSATE IN A TRAP 

In this section we consider a quasi one-dimensional 
condensate confined in an harmonic trap V cy ± = 



4 



tor: 




( Zr z 2 )/L 



FIG. 3: Energy shift in units of K 2 N m^L 2 z /h 3 V 2 for two im- 
purities in a homogeneous condensate with periodic boundary 
conditions. The interaction of BEC atoms is characterized by 
the dimensionless parameter K — gNo2mBL 2 /fi 2 V . The im- 
purities are located on the — 0.5-Lz, and 
z = 0.05L Z . 
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(25) 

(26) 

(27) 
(28) 



Here zb = \J h/mu>B is the ground-state width of the 1-D 
harmonic trap. By calculating the integrals in equation 
one finds 



A 



where 



2h ^ 



1 N exp 



i) r 



tt(z 2 + z b ) 2»v\ (29) 



x H v {zi)H v (z 2 ) 



z z b/( z B + z o) and h = z lN z B + z 0- Wc 

also have introduced the one-dimensional coupling con- 
stant kid = K/(27ra^_) with the radial confinement a\ = 
hjrriB'jJ^- The conditional level shift is shown in figurc[S] 
for different widths z§ of the impurity traps. As expected 
the shape of the curves coincides for distances larger than 
the ground state width of the impurity traps. Distances 
smaller than zq are excluded because we have assumed 
that there is no direct scattering interaction between the 
impurity atoms. 

It is interesting to note that different from the case of 
a condensate in a box, the force between the impurities is 
not always attractive. One recognizes that this is only the 
case if the distance is sufficiently small. If the distance is 
larger than a certain value, in our case i\ — z^ as 2 ■ 0.6, 
the force becomes repulsive. 

We now consider the case of a weakly interacting 1- 
D gas. In order to solve the Gross-Pitaevskii equation 
we make use of the Thomas-Fermi (TF) approximation. 
Although the results obtained in this way cannot be 
smoothly connected to the ideal-condensate case, the TF 
approximation allows to derive a compact expression for 
the level shift. The TF condensate wavefunction is given 
by 



FIG. 4: Influence of condensate geometry on energy shift in 
a box with periodic boundary conditions. A is in units of 
K 2 2ma/4li 3 , K is defined as in figure |3] The impurities are 
located on the z-axis with a distance of 0.5L Z . L z — 12 ■ 
10 -6 m, and L x = L„ has been varied. 



M*) = \— i- 

V fflD 



p2 



(30) 



where the TF radius is given by Rtf = y7/j/mBw|. \i 
denotes the chemical potential and the one dimensional 
interaction parameter (?id is defined analogous to kid- 
To solve the Bogoliubov-de Gennes equations analytically 
further approximtions are needed as discussed in [12| ■ We 
here take over the results for the functions /• obtained 



mBw|z 2 /2. We first consider the case of an ideal, i.e. 
noninteracting gas. In this case the solutions of the 
Gross-Pitaevskii equation and the Bogoliubov-de Gennes 
equations are just the solutions of the harmonic oscilla- 



2j + l 



2R 



TF 



2[i 
~E~ 



Ft 2 



P, 



R 



TF . 



(31) 
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FIG. 5: Energy shift A in an ideal 1-D condensate as function 
of impurity distance for different ratios of zb/zo- The shift is 
given in units of ljbt (z b + Zo) ■ The inset shows 

a magnification for small distances. Here z\ — — Z2- 



As shown in Appendix C the result of equation i|33|) 
should be valid as long as the following conditions are 
fulfilled 



i?TF 



i?TF 



3> max 



Z2 ^ /c 

>>n n 



~0 



Rtf J 



(35) 
(36) 



Here, Sr denotes the distance of one of the impurities 
to the edge of the condensate and ( = hajB/^/i is the 
Thomas-Fermi parameter. Furthermore the interaction 
strength of the condensate has to fulfill the condition 



9m > 



2Hlu^ Rtf 



3N 

Hence, we have the restriction 



A <C min < ojb 



3AWd 
I6h 2 R^pUjB 



(37) 



(38) 



with the energies Ej = Hujb \J j(j + l)/2. The Pj are Leg- 
endre polynomials. Using the completeness of the Legcn- 
dre polynomials one can explicitely evaluate expression 
(TTHj) . One finds that the sum including the J = term 
vanishes if the overlap of the impurity wavefunctions is 
negligible: 



E 



(32) 



2fflD 



dz \<f> (z - zi)| 2 \4> (z - z 2 )\ 2 ~ 0. 



Thus the energy shift l|18|) is determined only by the j = 
term which yields the simple expression : 



A = 



"ID 



8fti?TF5lD 



(33) 



This result does not depend on the distance of the impu- 
rity atoms which is due to the Thomas-Fermi approxima- 
tion. The shift is always positive and becomes larger for 
smaller interactions in the BEC and for larger ID con- 
finement. It is instructive to express A in terms of the 
impurity-BEC stattcring length a; and the BEC scatter- 
ing length a. One finds 



A = 



ujj_ m\ 



4 mf i?TFO 



(34) 



Here we have m\ = 2mB?ns / (mb + ms)- Thus assuming 
a tight transversal confinement with uj± = 2ir x 10 4 Hz, 
a large scattering length between impurities and BEC 
a; = 200nm, a small scattering within the BEC a = 5 
nm, a small trap with i?TF = 20^m, and ms ~ tob one 
finds a conditional frequency shift of 2-7r x 10 3 Hz. 



V. CONCLUSIONS 

In the present paper we have analyzed the interaction 
of impurity atoms in a Bosc-Einstein condensate local- 
ized at specific positions by tight confining potentials. It 
was shown that in addition to the level shift caused by 
s-wavc scattering with the macroscopic condensate field 
there are also contributions from the interaction with vac- 
uum fluctuations of the Bogoliubov phonons. The self- 
and conditonal energy shifts were calculated for a BEC in 
a box with periodic boundary conditions. It was shown 
that size and sign of the conditional energy shift depends 
on the separation of the impurities and is largest for a 
highly anisotropic condensate geometry and for small in- 
teractions within the condensate. With increasing inter- 
action of the condensate atoms the spatial dependence 
becomes less and less pronounced. Motivated by these 
findings the level shift in a quasi one-dimensional har- 
monic trap was calculated. In the Thomas-Fermi limit 
a rather simple analytic expression was obtained from a 
Bogoliubov approach. For small trap sizes a conditional 
frequency shift in the range of several kHz seems feasible 
which could be of interest for the implementation of a 
quantum phase gate. 
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APPENDIX A: DERIVATION OF THE 
EQUATION OF MOTION FOR THE 
STATISTICAL OPERATOR 

The total statistical operator of both the condensate 
and the impurities is denoted by %. Its time evolution 
is then given by the Liouville-von Neumann equation 



mm = H,x(t) 



, where H = H B + H$ + flint is 



the Hamiltonian of the whole system, with H B being the 
Hamiltonian of the condensate, H$ that of the impurities 
and H m t the interaction. Changing into the interaction 
picture yields 



ihd t x(t)= H int (t), X (t) 



Formal integration and resubstitution leads to 



(Al) 



mm = H int (t),x(t ) 



lh Jt Q 



H int (t), H- mt (t'),x(t') 



(A2) 



Here, to is the time when the interaction starts. The sta- 
tistical operator for the impurity atoms can be obtained 
by tracing out the condensate, i.e. g(t) = Ttb [jt(f)]- 
This yields 



iW t g(t) =Tr B H iat (t),x(t ) 



dt' Tr B 



to 



flint (*), H int (t'),x(t') 



Following the standard approach we assume that the in- 
fluence of the impurity atoms on the condensate can be 
neglected and that the statistical operator of the whole 
system seperates as 

X(t) = g{t) ® g B (t) + x COT r(0 « g(t) <8> g B (h) • (A4) 

Furthermore since we have incorporated the mean-field 
contribution to the free Hamiltonian of the impurities, 
the expectation value of the interaction Hamiltonian van- 
ishes. , i.e. Ti'e (^gB{to)H mt (t)j = 0. With these approx- 
imations we obtain 



d t ~Q{t) = - p jT dt' 

• Tr R ( [flirt (*), H int {t'), g{t') ® g B (t ) 



(A5) 

The interaction Hamiltonian in the interaction picture 
can be expressed as 



H int (t) = Y, t)(a0,t\ (y^iC*) + ^§B 2 {t) 



a, 



where 



(A6) 



xp,t){aP,t\ = e*(* 8+ * B )*|a/3) e-*(^ s +* B >. 



(A3) Substituting this into eq. (|A5|) yields 



(A7) 



dtQa0,i6(t) = J t eaf3,-f8{t') (t)Bt (t') ^ - K a K^) + ^fli(i)fl 2 (t')) ( K a K /3 ~ K/3« 7 ) 

+ ^fl 2 (t)fli(t')) (K a np - Ka K S ) + (B 2 (t)B 2 (t')} (4 - npn s ) 
+ (B x {i!)Bi{f)\ (tf - K aKl ) + ^Si(t')fl 2 (<)) {k 7 k s - Ki3 k 7 ) 

+ (fl 2 (t')fll (*)) ( K -( K S ~ K aKS ) + (jhfflBiitj) (4 ~ KpKs) 



(A8) 



r 



APPENDIX B: BOGOLIUBOV THEORY 



In this appendix we briefly summarize the main results 
of the Bogoliubov approach. We start with the hamilto- 
nian of the Bose gas in s-wave-scattering approximation 



The field operator tp of the condensate is then devided 
into a C-numbcr function ?po which represents the con- 
densed part of the Bose-gas and an operator £ of quantum 
fluctuations: ip(r) = tp (r) + £(r). The wavefunction of 
the condensate is given by the Gross-Pitacvskii equation 



fl B = J dr^t(r) 



2m B 



A + V ext (r)-n U(r) 



f / dr^(r)^(r)ip{r)^(r) 



(Bl) 



" A + V^ t (r)-n + g\Mr)\ 2 )Mr) = 0- (B2) 



2m B 
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By plugging this into the Hamiltonian and neglecting 
terms of the order 0(£ 3 ) and higher one gets 

H B ^H B +Jdr ||t( r ) (~^A + Kxt(r) - m) l(r) 

+ |(4|Vo(r)| 2 ^(r)|(r) 

+^ 2 (r)f(r)et(r) + ^o 2 (r)e(r)e(r) 



The terms linear in £ vanish because of the Gross- 
Pitaevskii equation. The term does not depend on 
operators and is without consequence. In order to diago- 
nalizc the Hamiltonian we employ the Bogoliubov ansatz 

i(r) = J2'Mr)b v -v:(r)bt (B4) 
i\v) = Y^nl{vp v -v v {T%. (B5) 

V 

Here, b\, and b v are bosonic creation and anihilation op- 
erators of the Bogoliubov quasi-particles. The prime at 
the sum indicates that the ground state is excluded in 
the summation. If the wave functions u v and v v fulfill 
the Bogoliubov-de Gennes equations (V>o is taken to be 
real) 



h 2 A 

h 2 A 



+ Kxt(r) - H 



Uu+g\ipo\ (2u v — v v ) = E v u v 

(B6) 



2?7lB 



+ Kxt(r) - fx 



v v + g\ipo\ (2v v -u v ) 



-E v v v , 



(B7) 



with the normalization 

{u v {v)ul,{v) - Wl/ (r)i£,(r)}dr = <W (B8) 
{v v {v)u v ,{T)-u v {r)v v .{r)}6.T = 0, (B9) 
the Hamiltonian takes the very simple form 
H B =H°-Y,' E„ I \v v {v)\ 2 Av+Y^ E v blb v . (BIO) 

With this the operators £ in the interaction picture can 
easily be calculated 



CM) = 2^ u v (r)b v e 



iE„t/h 



t +iE v t/h 



(Bll) 



APPENDIX C: VALIDITY OF EQ. l|33|> 

In order to estimate the range of validity of the ex- 
pression for the conditional shift in TF approximation, 



cq. (|33() , we start with the expression (see also eq. 1111) 



gj(l,2) 



M 

£ Ej 

j=0 J 



M 

E 



~n E 3 J-Rt 



dz \M*-*i)\ M*)f7(*) 



(CI) 



(B3) where /• = Uj — Vj. By using (|3 1|) we find 



A I 



3=0 



Rtf 



Rtf 



x \M*-*i)\ 2 \M*'-**)ffp ' 



Rtf Rtf 



(C2) 



5 where we have introduced 

M 



/ P M (x, x>) = ^ ^±ip„(x)P„(^) . (C3) 



n=0 



If Af — * oo the sum approaches the ^-function and we 
obtain equation i|32|) . On the other hand the solutions 
STfy of the Bogoliubov-de Genne equations used here are 
only valid for [l2l | 



<5r 
Rtf 



3> max 



v / m(m + i)C 



V2C 



v/2 'V y/M(M + 1) 



(C4) 



where <5r is the distance from the edge of the condensate. 




FIG. 6: Picture of /p for M = 20. 

This implies M < V^St/RtfC an d with 5r > .RtfVC 

also following from eq. (|C4|) we arrive at M <C y| . Thus 

the limit M — + oo cannot be taken in (|C3|) . Neverthe- 
less even for a finite but sufficiently large upper limit of 
summation M the sum is to a good approximation zero 
as can be seen as follows: In figure El/p° is shown. One 
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recognizes a pronounced central maximum. The first in- 
tegral over z in equation 1|C2|) only contributes if there is 
an overlapp of the maximum of f^ 1 and the ground state 
<f>o (z — zi). The same holds for the second integral over 
z' and <fio(z' — Z2). Hence, equation (|C2() vanishes if the 
distance of the impurities is much bigger than the width 
of the central maximum. We thus need to estimate the 
width of this central peak. With the Stirling formula one 
finds asymptotically for large (and even) M 

/ P M (0,0)«-. (C5) 



Since / /p^ (0, s) ds = 1 the width of the central peak can 
be approximated as As = tt/M. This finally yields the 
condition 

for which the sum in l|C2|) is approximately 0. It should 
be noted that we have assumed the Thomas-Fermi limit 
( < 1, which is essential for the analytic solution of the 
Gross-Pitaevskii and Bogoliubov-de Gcnnes equations. 
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